Lagrangian statistics in fully developed turbulence 
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The statistical properties of fluid particles transported by a fully developed turbulent flow are 
investigated by means of high resolution direct numerical simulations. Single trajectory statistics 
are investigated in a time range spanning more than three decades, from less than a tenth of the 
Kolmogorov timescale up to one large-eddy turnover time. Acceleration and velocity statistics show 
a neat quantitative agreement with recent experimental results. Trapping effects in vortex filaments 
give rise to enhanced small-scale intermittency on Lagrangian observables. 



PACS numbers: 47.27.-i, 47.10. +g 

The knowledge of the statistical properties of particle 
tracers advected by a fully developed turbulent flow is 
a key ingredient for the development of stochastic La- 
grangian models in such diverse contexts as turbulent 
combustion, pollutant dispersion, cloud formation and 
industrial mixing |jj,|2|,y,|4j]. Given the importance of this 
problem, there are comparatively few experimental stud- 
ies of turbulent Lagrangian dispersion. Progress in this 
direction has been hindered mainly by the presence of 
a wide range of dynamical timescales, an inherent prop- 
erty of fully developed turbulence. Indeed, in order to 
obtain an accurate description of the particle statistics it 
is necessary to follow their paths with very high resolu- 
tion, i.e., well below the Kolmogorov timescale r^, and 
for a long time lapse - of the order of an eddy turnover 
time Tl. The ratio of these timescales can be estimated 
as Tl/tti ~ R\ where the microscale Reynolds number 
R\ ranges in the hundreds for typical laboratory exper- 
iments. A recent breakthrough has been made by La 
Porta et al 0, Q who, borrowing techniques from high- 
energy physics, were able to track three-dimensional tra- 
jectories with a resolution of 0.05r^ and thus to study 
the statistics of particle acceleration. However, owing 
to the small measurement volume, trajectories could be 
followed only up to a few t v , preventing any investiga- 
tion of the long time correlations along particle paths. 
Conversely, the acoustical technique adopted by Mor- 
dant et al pj enabled them to successfully track parti- 
cles for durations comparable to Tl but could not ac- 
cess time delays of the order of t v , and was restricted to 
one-dimensional measurements. More conventional tech- 
niques, e.g. based on CCD cameras 0, reveal useful only 
for moderate Reynolds numbers (R\ ~ 100) due to their 
limitations in the acquisition rate. In addition, all labo- 
ratory techniques require a considerable amount of signal 
processing to get rid of experimental noise. Besides these 
drawbacks, increasing difficulties are met when attempt- 
ing to deal with multi-particle tracking. 



In the light of the above remarks Direct Numeri- 
cal Simulation of fully developed turbulence represents 
a valuable alternative tool for the investigation of La- 
grangian statistics [l[ll[lJ[ll[l!Qmm[l3and 
its important contribution has been recently reviewed in 
Ref. [4j. Nonetheless, computational approaches have to 
face three demanding requirements: (i) have the largest 
possible value of R\ in order to maintain the flow in a 
fully developed turbulent state; (ii) resolve properly the 
dissipative scales with the aim of computing accurately 
small-scale observables; (iii) follow particle trajectories 
for times comparable to a large-eddy-turnover time. High 
resolution and huge computational resources are there- 
fore needed to accomplish such a goal, making massive 
parallel computing by far and away the most appropriate 
tool. In this Letter we report the results of Direct Numer- 
ical Simulations of Lagrangian transport in homogeneous 
and isotropic turbulence on the parallel computer IBM- 
SP4 at CINECA, on 512 3 and 1024 3 cubic lattices with 
Reynolds numbers up to R\ ~ 280, a very accurate res- 
olution of dissipative scales, and a long integration time 
T « Tl- The Navier-Stokes equations 

dtu + u ■ Vtt = — Vp + i/Au , (1) 

are integrated on a triply periodic box by means of a fully 
dealiased pseudospectral code (the numerical parameters 
are listed in Table UJ). Energy is injected at an average 
rate e by keeping constant the total energy in each of 
the first two wavenumber shells 18]. With the present 
choice of parameters the dissipative range of lengthscales 
is extremely well resolved. With regard to the scaling 
properties of the velocity field, our results are in per- 
fect numerical agreement with previous numerical simu- 
lations at comparable R\ (see e.g. |l9]). Upon having 
reached statistically stationary conditions for the veloc- 
ity field, millions of Lagrangian tracers have been seeded 
into the flow and their trajectories integrated according 
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Rx Urms £ v y L T L t v T 5x N 3 Np 

183 1.5 0.886 0.00205 0.01 3.14 2.1 0.048 5 0.012 512 3 0.96 10 6 
284 1.7 0.81 0.00088 0.005 3.14 1.8 0.033 4.4 0.006 1024 3 1.92 -10 6 

TABLE I: Parameters of the numerical simulations. Microscale Reynolds number R\, root-mean-square velocity u rm s, energy 
dissipation e, viscosity v, Kolmogorov lengthscale 77 = (^ 3 /e) 1//4 , integral scale L, large-eddy turnover time Tl = L/u r ms, 
Kolmogorov timescale r v = (ly/e) 1 ^ 2 , total integration time T, grid spacing Sx, resolution iV 3 , and the number of Lagrangian 
tracers N p . 




FIG. 1: Snapshots of a subset of Lagrangian particles ( m 12000) at three different times t = (0.2, 0.5, 0.8)Tt. Gray scales 
encode different velocity amplitudes. 



to 



— =u(X(t),t), 



(2) 



over a time lapse of the order of Tl (See Fig. Parti- 
cle positions X(t) and velocities v(t) = X(t) have been 
stored at a sampling rate 0.07T,,. The forces acting on the 
particle - pressure gradients Vp(X(t), t) , viscous forces 
vAu(X(t), t) and external input - and the resulting par- 
ticle acceleration a(t) = v(t) have been recorded along 
the particle paths every 0.14r^. The resulting database 
permits detailed study of the statistics of Lagrangian ve- 
locity and acceleration over a range of timescales span- 
ning more than three decades. 

In what follows, we will focus on the description of the 
statistics of single particle trajectories, deferring the dis- 
cussion of multi-particle statistics to a forthcoming pub- 
lication. The analysis reveals the presence of frequent en- 
trapment events within vortical structures (see Fig. [21 . 
The characteristic frequency of these spiraling paths is 
comparable to r" 1 and their duration can be as long 
as 10t v . The velocity experienced during these events 
can attain 5M rms with an ensuing acceleration 5u rms /T n 
as large as 80a rms . As we will show below, these events 
have a significant impact on time correlations up to 1Qt v . 

In Fig. we present the probability density functions 
(PDFs) of velocity differences S T v = v(t + r) — v(t) taken 
along a given direction and rescaled by their variance 
((Stv) 2 ) 1 / 2 . For the shortest time delays shown in the 
figure the velocity increment PDF approaches the distri- 
bution of the acceleration a(t) which is characterized by 
a strong intermittency, stretched exponential tails and a 
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FIG. 2: Trajectory and time series. Left panel: 3D trajectory 
of a trapping event in vortex filament. Acceleration and ve- 
locity fluctuations here reach about 30 and 2 r.m.s. values, 
respectively (right panels). 



flatness (a 4 ) /(a 2 ) 2 « 40. Classical Kolmogorov scaling 
arguments (see e.g. Ref. |2(j) yield for the acceleration 
variance the prediction 



{a 2 )=a s^v- 1 ' 2 



(3) 



We measure the constant an = 3.5 ± 0.3, in agreement 
with the results of Ref. [13| at comparable R\. Yet, 
one needs to notice that intermittent corrections may 
introduce in the definition of eto some dependency on 
Reynolds number 0. At larger time separations the 
PDFs are decreasingly intermittent and eventually be- 
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FIG. 3: Probability density functions of velocity increments 
and acceleration, normalized with their variance. Curves refer 
to time increments r = (97, 25, 6, 0.7)t v from inside to out- 
side, and to the acceleration (outermost curve). In the inset, 
the kurtosis K(t) = ((Si-?; 4 ) /{{5 t v) 2 ^ 2 for the entire time in- 
terval 0.07t v -t- 2Tl- The saturation of K(t) at small time 
increments is an indication of the high numerical resolution. 



come slightly sub-Gaussian for r w Tl, with a flatness 
((d T v) A ) I {(S T v) 2 ) 2 2.8. Assuming a scaling law for 
the second-order Lagrangian velocity structure function 
in the time-range t v -C t -C Tl, dimensional analysis 
|2(ij . predicts {{S T v) 2 } — CqEt. Remark that this predic- 
tion is not affected by intermittency corrections owing 
to its linear dependence on the energy dissipation rate 
[2ll E2I E3 | . Therefore, Co is expected to be universal. 
For our data the scaling behavior holds on a relatively 
narrow window in r, leading to some uncertainty in the 
determination of the numerical prefactor. We estimate 
Co = 5 -7- 6, not far from previous measurements 

The intermittency in the PDFs of 5 T v can be conve- 
niently quantified in terms of the Lagrangian velocity 
structure functions S p (t) = ((S T v) p ), that are expected 
to behave as power laws r^ p . The scaling exponents 
£ p can be evaluated by looking at the logarithmic slope 
d log Sp(r)/d log r that should display a plateau in the 
range r ?) <C r <C Tl. From Fig.01it is clear that it is very 
difficult to extract the values for £ p . However, by means 
of the extended self-similarity procedure [25j, it is possi- 
ble to estimate the relative exponents ^4/^2 = 1.7 ±0.05, 
£ 5 /£ 2 = 2.0 ± 0.05, £ 6 /£ 2 = 2.2 ± 0.07, in fair agreement 
with those obtained in Ref . Q ■ The range of time de- 
lays over which relative scaling occurs is IOt^ to 70r r) . In 
this range anisotropic contributions induced by the large 
scale flow appear to influence the scaling properties. 

It is interesting to remark that for values of r in the 
range from t v to 10t v , the local slopes are significantly 
smaller and tend to accumulate around the value 2 for all 
orders. This relevant correction to scaling cannot be at- 
tributed to the influence of the dissipative range r < t v , 



since the latter would increase the value of the local slope, 
rather than decreasing it. A similar effect can be detected 
in Eulerian structure functions as well, yet the intensity 
in the latter case is much less pronounced. These strong 
deviations in the Lagrangian scaling laws are most likely 
due to the trapping events depicted in Fig. |21 Indeed, 
the long residence time within small-scale vortical struc- 
tures introduces an additional weighting factor that en- 
hances the effect with respect to Eulerian measurements. 
Further insight into the single particle dynamics can be 
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FIG. 4: Log-log plot of Lagrangian structure functions of or- 
ders p = 2, 4, 6 (bottom to top) vs r. Bottom right: logarith- 
mic local slopes d log S p (r) / d log r (same line styles) . Top left: 
relative local slopes with respect to the second order struc- 
ture function dlog S p (r)/d log £?2(t), for p = 4, 6. Data refer 
to the v x component. The two other velocity components 
exhibit slightly worse scaling due to anisotropy effects. 

gained by investigating the various force terms concur- 
ring to determine the total particle acceleration 



d 2 X{t) 
dt 2 



= -V P (X, t) + vAu(X(t),t) + F(X(t), t) 



As shown in the inset of Fig. [SI it is evident that the 
main source of acceleration is the pressure gradient 01 j 
whereas the viscous term is important only near rela- 
tively strong acceleration events and the forcing contri- 
bution is negligible. A more quantitative test is given 
by the comparison of the PDFs of the three forces in 
Fig. [3J To conclude, we have presented the analysis of 
single-particle statistics in high Reynolds number flows. 
Our results fit well with previous experimental measure- 
ments. At variance with experiments, we can investigate 
the statistical properties of millions of particles on a wide 
range of time intervals, from a small fraction of the Kol- 
mogorov time up to the integral correlation time. We 
found clear indications that velocity fluctuations along 
Lagrangian trajectories are affected by multiple-time dy- 
namics. In the interval 10r, ; < r < Tl we observed 
anomalous scaling for Lagrangian velocity structure func- 
tions. For frequencies of the order of r^ 1 we noticed 
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cussions and B. Devenish for a careful reading of the 
manuscript. 
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FIG. 5: Log-log plot of PDFs for —d x p, vAu x , F x . The exter- 
nal forcing is virtually negligible, and the main contribution 
to large accelerations is made by pressure gradients. Inset: 
a typical evolution of the three terms along a particular tra- 
jectory. The strongest signal is d x p (dashed line), while the 
viscous force is activated only as a subleading response to 
pressure gradients (solid line). The force contribution is in- 
distinguishable from zero. 



that velocity fluctuations are affected by events where 
particles are trapped in vortex filaments. Events with 
trapping times much longer than expected on the basis 
of simple dimensional analysis appear frequently. The 
main novelty of Lagrangian single-particle velocity fluc- 
tuations with respect to Eulerian fluctuations is prob- 
ably the capture of trajectories induced by small-scale 
structures. For example, the event analyzed in Fig. El 
would have a much smaller weight in an Eulerian analy- 
sis because of large-scale sweeping past the fixed probe. 
Among the most challenging open problems arising from 
our analysis are how to incorporate such dynamical pro- 
cesses in stochastic modelization of particle diffusion Q 
and in the Lagrangian multifractal description |22l Eij ] . 

The simulations were performed within the keypro- 
ject "Lagrangian Turbulence" on the IBM-SP4 of Cineca 
(Bologna, Italy). We are grateful to C. Cavazzoni and 
G. Erbacci for resource allocation and precious techni- 
cal assistance. We acknowledge support from EU un- 
der the contracts HPRN-CT-2002-00300 and HPRN-CT- 
2000-0162. We also thank E. Leveque for useful dis- 
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